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ABSTRACT 



Aims. We look for a simple analytic model to distinguish between stellar clusters undergoing a bimodal hydrodynamic solution from 
those able to drive only a stationary wind. Clusters in the bimodal regime undergo strong radiative cooling within their densest inner 
regions, which results in the accumulation of the matter injected by supemovae and stellar winds and eventually in the formation of 
further stellar generations, while their outer regions sustain a stationary wind. 

Methods. The analytic formulae are derived from the basic hydrodynamic equations. Our main assumption, that the density at the 
star cluster surface scales almost linearly with that at the stagnation radius, is based on results from semi-analytic and full numerical 
calculations. 

Results. The analytic formulation allows for the determination of the threshold mechanical luminosity that separates clusters evolving 
in either of the two solutions. It is possible to fix the stagnation radius by simple analytic expressions and thus to determine the fractions 
of the deposited matter that clusters evolving in the bimodal regime blow out as a wind or recycle into further stellar generations. 
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1. Introduction 

The discovery of young massive stellar clusters or super stel- 
lar clusters (hereafter SSCs) with typical masses between sev- 
eral times 10^ Mq to several times 10^ Mq and radii in the 
range 1 - 10 pc has unveiled the dominant mode of star for- 
mation in starburst and interacting ga laxies (see, for example, 
iMcCradv et"!!] 120031; IWhitmorell2006 and references therein). 
SSCs dominate the structure of the ISM in their host galaxies 
via their large UV photon output and their high velocity out- 
flows or star cluster winds. Furthermore, the interaction of the 
winds from nearby SSC may lead to the inner structure of the 
large-scale galactic outflows, or super galactic winds, which con- 

, nect starburst events with the intergalactic medium as in M82 

. jTenorio-Tagle et al.ll2003[) . 

Here, we center our attention on the hydrodynamics of 
the matter reinserted within the SSC volume, particularly in 
the case of massive and compact SSCs. The latter have been 
shown t o be exposed to the strong radiative cooling (see 
ISilich et al. 2004) th at lead s to a bimodal hydrodynamic solution 
dTen orio-Tag le et al] l2007h . This phenomenon occurs because 
the cluster mechanical luminosity and the injected gas density 
are nearly linear functions of the cluster mass, while the cooling 
rate is proportional to the square of the density, and therefore if 
one considers very massive and c ompact clusters, ra diative cool- 
ing would always dominate (see ISilich et al.l 12003 for details). 
Radiative cooling depletes the energy in the densest inner re- 
gions enclosed within the so-called stagnation radius. The gas 
velocity at the stagnation radius is zero km s Thus the mat- 
ter reinserted through stellar winds and supernovae within the 
central volume is accumulated to become thermally and/or grav- 
itationally unstable leading to further episodes of star formation 
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(iTenorio-Tagle et ani2005h . On the other hand, the matter in- 
jected in the outer layers, between the stagnation radius and the 
cluster surface, composes a stationary wind. In the adiabatic so- 
lution and for low mass and extended clusters, the stagnation ra- 
dius is always found at the star cluster center However, if radia- 
tive cooling is taken into consideration, it progressively moves 
towards the star cluster surface as one considers more massive 
and compact clusters. 

In this work, as in IChevalier & Cleg3 (1198 5'), the density, 
temperature and velocity are functions of the radial coordinate 
only. Such an approach does not take into account the effects 
of the turbulence generated by individual stellar winds and su- 
pernovae within the cluster volume. By changing the density, 
velocity and temperature distributions, it is thought to provide 
an additional transport of energy towards the smallest condensa- 
tions. However it cannot prevent catastrophic cooling within the 
densest parts of the flow because it results from the net energy 
balance inside the cluster 

The properties of such clusters are described here by an ap- 
proximate analytical formulation, in agreement with our former 
semi-analytic and numerical calculations (see sections 2 and 3). 
We display our results in a series of universal diagrams valid for 
all massive and compact SSCs. Section 4 deals with the impact 
of different heating efficiencies. 

2. Adiabatic vs radiative solutions 

The original framework to study the fate of the matter deposited 
by supernovae and stro ng stellar winds within a star cluster vol- 
ume was developed bv iChevalier & Clegd (11985 ). In their adia- 
batic model the kinetic energy released by individual sources is 
completely thermalized in situ via direct collisions of the high 
velocity gaseous streams from neighboring sources. This results 
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in a high central temperature (~ 10^ - lO^K) and pressure that 
exceeds the pressure of the surrounding interstellar medium. The 
velocity of the injected matter then acquires a particular distri- 
bution, increasing from zero km s"' at the star cluster center to 
the speed of sound (csc) at the star cluster surface, approaching 
rapidly a supersonic value (Vaoo) as it moves outside the cluster 
volume. Csc is defined as: csc = [(y - + ^)V^^Vaco', where 
Vaco = [2Lsc/Msc]^^^ is the adiabatic wind terminal velocity, 
Lsc is the star cluster mechanical luminosity, and Msc is the 
rate of mass injection provided by supernovae and stellar winds. 
The condition for a stationary flow is then: 



Msc - ^^^IcPscCsc, 



(1) 



where Rsc is the radius of the cluster and psc is the density of the 
outflow at the star cluster surface. The relation holds whenever 
Psc reaches the proper value, which is in direct proportion to the 
density at the stagnation radius (i.e. in the center), p^,: psc - 
[(r + l)/(6r + 2)](''''+i)/(5r+i)p^^^ ^here p,, is (ICanto et al.(l2000l) 
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Here q,„ - (3Msc)/(^^Rlc^ mass deposition rate per unit 
volume which for a given Rsc scales linearly with the star cluster 
mass, Msc, and y is the ratio of specific heats. Note that in the 
adiabatic solution, if one fixes the radius of the cluster, Rsc, and 
the adiabatic terminal velocity, Vaco, the central density and the 
density at the star cluster surface grow linearly with the mass of 
the considered cluster, since the stagnation radius is located at 
the star cluster center and thus the density profiles remain self- 
similar. 

It has recently been shown, however, that the adiabatic model 
is inadequate for very massive and compact star clusters, as 
in these cases ra diative cooling becomes a dominant factor 
dSilich et alj|2004l) . The authors demonstrated that in the radia- 
tive case the wind density at the stagnation radius changes to 
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where jt/,- = 14/1 Inifj is the mean mass per ion (niH is the proton 
mass), A{Tst,Z) is the cooling function (the cooling rate is Q - 
ngUiMTst, Z)), Tsi and c,, are the temperature and speed of sound 
at the stagnation radius, and Z is the metallicity of the plasma. 
The parameter rf describes how much of the deposited energy 
is transformed into energy of the outflow (see §4 for a detailed 
discussion). We assume rj - I everywhere in §2 and §3. 

For clusters of a given size, the density at the stagnation ra- 
dius grows in all cases with the mass of the cluster. In the case of 
low mass clusters it results from a combination of two effects: i) 
the mass deposition rate per unit volume, <7,„, is a linear function 
of Msci~ Msc)- ii) the temperature at the stagnation radius, Tj,, 
decreases with increasing cluster mass Msc making Cj, smaller 
and thus the term in brackets in the rhs of Eq. ([3]) larger. The 
combined action of these two factors then causes that the radia- 
tive solution proposed by Silich et al. (2004) is almost identical 
to the adiabatic solution of lChevalier & CleggI (Il985b . 

However, if one considers more massive clusters (or clusters 
with a higher mechanical luminosity, Lsc) one would soon reach 
a critical value for which strong radiative cooling would set in at 
the star c luster center, drastically c hanging the inner wind struc- 
ture (see lTenorio-Tagle et alj2007l) . For clusters above this limit, 
the temperature r„ at the stagnation radius Rg, does not change 



any more and remains equal to the value that corresponds to the 
highest possible pressure there (see section [3T| below). Thus the 
second term in Eq. (O remains fixed and the density at the stag- 
nation radius, p,,, grows as the square root of the star cluster 
mass, whereas the mass deposition rate continues to grow lin- 
early with the mass of the considered cluster This implies that 
above the threshold line, the stationary condition Eq. ([1]) cannot 
be fulfilled unl ess the stagnation radius mo ves towards the star 
cluster surface. iTenorio-Tagle et al] ( l2007l) showed how to ex- 
tend the semi-analytic approach proposed in lSilich et alj (|2004|) 
for this case. The authors also found the semi-analytic solution 
for the outflow beyond the stagnation radius, and compared it 
with full numerical calculations performed with the finite dif- 
ference hydrodynamic code ZEUS for which the cooling rou- 
tine was modified to make it suitable for the modeling of ex- 
tremely fast cooling regions. In particular, they include the cool- 
ing rate in the computation of the time-step assuming that the 
amount of energy that can be radiated away from a given cell 
during one time-step must be smaller than 10% of its internal 
energy. The time-step is decreased to meet this cooling rate con- 
dition. However, since this could lead to extremely small time- 
steps, which would substantially degrade the overall code per- 
formance, they do not allow the global time-step to decrease 
below 0.1 times the "hydrodynamic" time-step determined by 
the Courant-Friedrich-Levi criterion. If a certain cell requires 
an even smaller time-step due to the cooling rate condition, the 
time-step was subdivided even further and then was used to 
numerically integrate the energy equation only in the affected 
cell(s). This "time refinement" was applied only locally, so the 
CPU time is not wasted in cells where the high time resolution 
is not required (see .Tenorio-Tagle et al...2007. for more details). 



3. The catastrophic cooling regime 

3.1. The temperature at the stagnation radius 

Eq. ^ allows one to determine the pressure at the stagnation 
radius f - pstkTsi/fi where k is the Bolt zmann constant, and 
U — 14 /23OTtf is the mean mass per particle. pTenorio-Tagle et all 
(I2007h showed that in the bimodal regime the pressure at the 
stagnation radius always attains the maximum value allowed by 
the shape of the cooling function, and thus one can find the cor- 
responding temperature numerically if Vam and the shape of the 
cooling function, A(r, Z), are known. 

Fig.[T]displays Tj, for different values of the adiabatic wind 
terminal speed (or different ratios of energy to mass deposition 
rates) and different metallicities of the thermalized plasma. The 
temperature at the stagnation radius grows approximately as a 
quadratic function of Vaoo and remains always below the adi- 
abatic value (dotted line in Fig. [1). For practical purposes one 
can suggest that it does not depend on the plasma abundance 
(see Fig. [TJ. Note that f or these calculations and throughout the 
paper we have used the Ray mond. Cox. & Smith ' (1916) equilib- 
rium cooling function tabulated and updated bv iPlewa (.1995.) . 



3.2. The threshold mechanical luminosity 

In the case of a homogeneous stellar mass distribution, the equa- 
tion of mass conservation 
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Fig. 1. The temperature at the stagnation radius as a function 
of the adiabatic wind terminal speed, Vaco = 2Lscl Msc- Solid 
and dashed lines represent the temperature at the stagnation ra- 
dius for solar and 10 times solar metallicities, respectively. The 
adiabatic wind stagnation (central) temperature is shown by the 
dotted hne. 
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Fig. 2. The ratio of the density at the star cluster surface to 
that at the stagnation radius in the catastrophic cooling regime. 
Solid, dotte d and dashed lines present semi-analytic results (see 



ISilich et al.l l2004) for Vaco = 1000km s 
and Vaoo — 2000km s"', respectively. 
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can be easily integrated. At the star cluster edge this results in: 
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As mentioned above, in the adiabatic case the density at the 
star cluster surface scales linearly with that at the central stagna- 
tion radius. This is not exactly the case in the radiative solution. 
Nevertheless, psc/Psi remains restricted into a very narrow in- 
terval, 0.28 < psclPsi ^ 0.42, (see Fig.|2]l, which allows one to 
propose a simple analytic expression for the stagnation radius in 
reasonable agreement with the results found in semi-analytic and 
full numerical calculations. Indeed, assuming by analogy with 



the adiabatic solution, that psc - apst, where a is a fiducial co- 
efficient, one can combine Eqs. ^ and (O to derive: 
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In Eqs. (|6]l and ^ A,, and c,, are the cooling function value 
and speed of sound at the stagnation radius, respectively. To de- 
rive these equations we have used the adiabatic relation for the 
sound speed at the star cluster surface: c^^ - {y- 1)^aoo/(')'+ !)■ 
When the mechanical luminosity of the cluster exceeds Lent, the 
stagnation radius detaches from the star cluster center, which im- 
plies that Lcrii is the threshold luminosity which separates clus- 
ters evolving in the catastrophic cooling (bimodal) regime from 
those evolving either in the quasi-adiabatic or in the radiative 
regimes, with their stagnation radius located at the star cluster 
center. Note that for the latter cases Eq. (|6]l is not valid. 



3.3. The analytic formulation 

The analytic values for the threshold luminosity and the stag- 
nation radius have been compared with our semi-analytic re- 
sults and full numerical calculations (see Fig. |2] and [3] in 
iTenorio-Tagle et al]|2007l) . For that we have used first Fig. [T]to 
obtain the temperature at the stagnation radius for a given value 
of the adiabatic wind terminal speed, Vaoo. We then use this tem- 
perature to calculate the sound speed, c,,, and the corresponding 
cooling function value. A,,, and then use these values to cal- 
culate the threshold luminosity and the stagnation radius from 
Eqs. (|7|i and (|6]l. The agreement between the approximate ana- 
lytic expressions and semi-analytic results is good if the fiducial 
coefficient a - 0.28 (see Fig.[3]l. 

Using Eq. (|6]l one can also estimate the fraction of the de- 
posited matter that clusters evolving in the bimodal regime re- 
turn to the ISM of their host galaxy: 



Ri 



Ms 



Ri 



Lcril 



1/2 



(8) 



This quantity decreases monotonically with the star cluster 
mechanical luminosity (see Fig. O as wa s also found in the 
ID hy drodynamic simulations performed bv lTenorio-Tagle et al.l 
(I2007h . 

The analytic expression ([8]) reproduces the numerical results 
with an accuracy better than ~ 25%. Minor differences between 
the analytic and the numerical outputs indicate that psc/Psi de- 
viates from the assumed constant value. 



4. Heating efficiency 

Here we follow the IChevalier & Cl eg^ ('1985') formulation 
sup ported by full numerical calcu lations (see, for example 
iCanto et aljr2000t [Rockefeller et al.ll2004i) and assume that the 
injected matter is thermalized in situ via random encounters of 
high velocity stellar winds and supernovae ejecta. The efficiency 
of this process is an important, but poorly known parameter 
(Stevens & Hartwell 2003). This para meter has been discussed 
by different groups (see, for example , iBradamanteet all 1 19981 : 
iMelioli & de Gouveia Dal Pino! |2004| and references therein) 
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Fig. 3. The comparison of analytic results with semi-analytic and 
full numerical calculations. Panel a) displays the critical lumi- 
nosity calculated for three different values of the adiabatic wind 
terminal speed: Vaoo - 1000 km s"', Vaoo = 1500 km s"' and 
Vaoo - 2000 km s"'. Panel b) displays how the stagnation ra- 
dius shifts from the center when the mechanical luminosity ex- 
ceeds the critical value. Panel c) shows the fraction of the matter 
supplied by supernovae and stellar winds that the cluster returns 
to the ISM of the host galaxy. In aU panels solid lines present 
the analytic results. The results from the semi-analytic and full 
numerical calculations are shown by the dotted line and by the 
cross symbols, respectively. The fiducial coefficient a was set to 
a = 0.28 in all calculations. 

who argued for very different values of heating efficiency that 
vary from 100% to a few per cent. 

Note that different authors have different prescriptions for 
the physical processes that may lead to an incomplete transfor- 
mation of the deposited kinetic energy. In our view, the high 



density of the sources may cause an important fraction of the 
star cluster mechanical luminosity to be immediately radiated 
away. This would happen locally at the sites where the ejecta 
from nearby stellar winds and supernovae would interact, in- 
hibiting in this way the possibility of having such energy evenly 
spread throughout the cluster volume. Turbulence may also ab- 
sorb some fraction of the deposited energy. The immediate con- 
sequence of such a sudden depletion of energy is to reduce the 
thermal pressure gradient responsible for the acceleration of the 
injected matter and thus the magnitude of the resultant high ve- 
locity outflow. We have incorporated all uncertainties related to 
this process into a parameter rj that represents the fraction of 
the deposited mechanical energy per unit time which is finally 
evenly spread throughout the cluster volume and participates in 
the thermal acceleration of the injected matter. The heating effi- 
ciency then effectively reduces the value of the deposited energy, 
Lsc, but does not affect the mass deposition rate, Msc- 

Fig.|4]shows the impact of a reduced heating efficiency on the 
plasma temperature and how it drastically decreases the thresh- 
old luminosity. What the most reasonable value of the heating 
efficiency in super star clusters is, remains an open question in 
need of both detailed full hydrodynamic al calculations and care- 
ful analysis of the available observational data. Nevertheless, the 
threshold mechanical luminosities calculated for a reasonable 
range of heating efficiencies, 0.1 < tj < 1 (shown in Fig. ^p), 
suggest that many SSCs may be found well above t he threshold 
line, as in the case of the M82-A1 supercluster. See lSmith et al.l 
(12006 ). who, based on the size of its associated HII region, es- 
timated a heating efficiency rj < 0. 1 . A critical revision of the 
available super star cluster data will be the subject of a forth- 
coming communication. 

5. Conclusions 

We have provided a simple analytic formulation that allows for 
an easy determination in the energy deposition rate vs cluster 
size diagram, of the location of the threshold mechanical lumi- 
nosity that separates clusters evolving in either of the two pos- 
sible hydrodynamic solutions: the normally assumed negative 
feedback solution in which all mass supplied by SNe and indi- 
vidual stellar winds leaves the cluster as a cluster wind - 0) 
and the bimodal solution (0 < /?„ < Rsc) in which the reinserted 
matter accumulates within the volume defined by R^t- Simple 
analytic expressions, accounting also for a reduced heating effi- 
ciency, allow one to calculate the position of the stagnation ra- 
dius and the fraction of the deposited matter that clusters evolv- 
ing in the bimodal regime return to the ISM of their host galaxy. 
Our results are in good agreement with our former semi-analytic 
and numerical results. 
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